AD-A049  772 


UNCLASflFXCO 


PRINCETON  UNIV  N J •«$  DYNAMICS  LAB  m/%  %q/h 

^ fin-inouceo  Shock  wave  turbulent  b— etc(u) 
JUL  77  0 S OOLLXNBt  C 0 COSAO*  S M BOBOONOPF  OAABt9-76-BoOtB9 

ARO*i90t«.i-CX  NL 


3 -70 


i 


ABSTRACT 

Surface  pressure  distributions  and  oil  streak  patterns 
have  been  measured  in  the  three-dimensional  flowfield  caused  by 
the  interaction  of  a skewed  shock  wave  with  a two-dimensional 
turbulent  boundary  layer.  The  boundary  layer  was  developing  on 
a 45  sweptback,  sharp  leading  edged  plate.  Tests  were  made  at 

A 

a freestream  Mach  number  of  2.95,  a freestream  unit  Reynolds  number 

of  -6.3  X lO^/m'^  (1.6  X 10^^  in"*^),  a stagnation  temperature  of  about 

265^K  and  at  near  adiabatic  wall  temperature. 

Measurements  have  been  made  for  generator  angles  in  the 
^ ^ ^ 

range  0 to  I 2®,^  i n 2J^^  i ncrements  . These  measurements  represent 
the  first  phase  of  a study  which  will  include  heat  transfer  mea- 
surements and  details  of  the  flowfield.  The  measured  data  are 
presented  and  compared  with  data  obtained  under  the  same  shock 
wave  conditions  but  with  a different  boundary  layer. 
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NOMENCLATURE 


P local  static  pressure 

P calculated  static  pressure  ahead  of  oblique 

shock  wave 

P^  calculated  static  pressure  behind  oblique 

shock  wave 

X coordinate  parallel  to  tunnel  axis  measured  from 

the  shock  generator  leading  edge 

Y coordinate  normal  to  the  X axis  in  the  piane  of 

the  test  surface  measured  from  the  shock  generator 
leading  edge . 

X X coordinate  measured  from  the  calculated  shock 

wave  location  (see  Fig.  4). 

Oq  shock  generator  angle 

6 boundary 

boundary 
rise  for 


iayer  thickness 

layer  thickness  at  beginning  of  pressure 
given  value  of  Y . 
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I NTROOUCT I ON 


The  experimental  study  reported  briefly  on  herein  represents 
the  first  phase  of  a detailed  investigation  being  made  of  the  three- 
dimensional  flowfield  of  a v i scous- i nv  i sc i d interaction  between  a 
shock  wave,  from  a sharp,  45“  swept  back  leading  edged  generator, 
and  a turbulent  boundary  layer.  This  study  is  being  carried  out 
at  a Mach  number  of  2.95  and  a freestream  unit  Reynolds  number  of 
6.3  X 10^  m * (1.6  X 10^  in  *).  Data  obtained  from  this  study  will 
aid  in  the  modelling  of  this  complex  three-dimensional  phenomenon. 

The  current  program  is  an  extension  of  earlier  work  carried 
out  under  contract  to  the  Flight  Dynamics  Laboratory.  Detailed 
studies  were  made  using  the  thick  tunnel  wall  boundary  layer  (i.e., 

6 = .55")  at  a Mach  number  of  2.95  with  a sharp  leading  edged  gen- 
erator. For  this  case  the  Interaction  could  be  studied  up  to  about 
8 boundary  layer  thicknesses  away  from  the  shock  generator.  To 
obtain  a considerably  "larger  picture"  in  terms  of  6 a further 
set  of  tests  was  made  using  the  boundary  layer  which  developed  on 
a flat  plate  with  a sharp  straight  leading  edge.  Both  of  these 
test  programs  have  been  documented  in  detail  by  Oskam  et  at  (1975) 
and  Oskam  (1976). 

This  second  study,  using  a flat  plate  with  an  unswept  leading 
edge,  produced  a turbulent  boundary  layer  which  varied  in  thickness 
from  0.1  I"  to  0.18"  along  the  leading  edge  of  the  Interaction.  In 
terms  of  an  average  boundary  layer  thickness  however,  the  Inter- 
action could  be  studied  at  up  to  about  30  boundary  layer  thicknesses 
away  from  the  generator. 
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To  reduce  the  variation  of  boundary  iayer  thickness  along 
the  leading  edge  of  the  interaction,  an  important  factor  in  scaling, 
the  unswept  leading  edge  of  the  originai  study  has  been  replaced 
with  a 45°  swept  back  ioading  edge.  With  this  geometrical  arrange- 
ment, the  surface  flow  patterns  and  surface  pressure  distributions 
have  been  measured  for  generator  angles  from  2°  to  12°  in  2°  incre- 
ments. This  study  constitutes  the  first  phase  of  the  present 
contract . 

It  was  anticipated  that,  with  a swept  leading  edge  plate, 
the  boundary  layer  thickness  ahead  of  and  along  the  interaction 
direction  would  vary  considerably  less  than  for  the  straight  leading 
edge  model.  A comparison  of  the  two  sets  of  results  would  aid  in 
the  evaluation  of  the  importance  of  the  boundary  layer  thickness  as 
a scaling  parameter. 

A brief  description  of  the  experimental  facility  and  the  model 
is  given  in  Section  2.  The  measured  pressure  distributions,  isobar 
patterns  and  surface  flow  angles  are  presented  and  briefly  discussed 
in  Section  3.  Comparisons  of  the  swept  and  unswept  results  are 
made  in  Section  4.  This  brief  interim  report  provides  an  outline 
of  the  work  done  to  date  under  the  present  contract. 
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2.  EXPERIMENTAL  EQUIPMENT 

2 . 1 Wind  tunnel 

The  experimental  study  was  carried  out  in  the  Princeton 
Universiry  high  Reynolds  number  supersonic  tunnel.  This  tunnel 
has  a working  section  20cm  x 20cm  (8"  x 8"),  a nominal  freestream 
Mach  number  of  3 and  may  be  operated  at  stagnation  pressures  in 
the  range  4 x 10^  to  34  x 10^  Nm  ^ (60  to  500  psi).  Since  the 
compressed  air  supply  for  the  tunnel  is  stored  outdoors  the  fiow 
total  temperature  is  variable  but  is  typically  in  the  range  260- 
270®K.  The  actual  value  at  any  given  time  is  determined  by  the 
local  external  temperature,  the  J ou 1 e-Thomp son  drop  across  the 
main  regulator  valve,  and  the  amount  of  heat  transferred  from  the 
inlet  piping. 

In  the  current  study  al I of  the  tests  were  made  at  a stagna- 

5 -2 

tion  pressure  of  6.8  x 10  Nm  (100  psi)  corresponding  to  a free- 
stream unit  Reynolds  number  of  6.3  x 10^  m * (1.6  x 10^  in  *). 

With  the  above  mentioned  temperature  considerations  the  model  was 
at  near  adiabatic  wall  conditions  for  all  tests.  Under  the  above 

conditions  and  with  the  compressed  air  supply  tanks  at  a maximum 

6 “2 

pressure  of  approximately  20  x 10  Nm  (3000  psi)  the  tunnel  can 
be  run  for  about  6 minutes.  A photograph  of  the  tunnel  showing 
the  nozzle  and  test  section  is  given  in  Fig.  I. 

2 . 2 Model  configurations 

The  test  surface  consisted  of  a 1.27cm  (.5”)  thick  brass 
flat  plate  with  the  sharp  leading  edge  swept  back  45®.  This  plate. 
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approximately  67cm  (24")  long,  spanned  the  tunnel  horizontally  anu 
was  mounted  5.08cm  (2")  above  the  tunnel  floor.  To  avoid  any  pos- 
sibility of  flow  blockage,  the  lower  tunnel  wall  was  recessed  beneath 
the  plate. 

Mounted  veri-ically  between  the  test  plate  and  the  tunnel 
ceiling  was  the  1.27cm  (.5")  thick  sharp  leading  edge  shock  gener- 
ator of  length  25,4cm  (10").  To  avoid  leaks,  the  edges  adjacent  to 
the  test  plate  and  the  tunnel  ceiling  were  lined  with  nylon  to  ensure 
good  sealing  at  any  generator  angle. 

The  generator  is  supported  by  a streamlined  arm  projecting 
through  the  tunnel  side  wall  which  is  linked  to  a drive  wheel  through 
a calibrated  screw  system.  A line  drawing  of  this  arrangement  is 
shown  in  Fig.  2.  This  allows  the  generator  to  be  manually  adjusted 
to  any  chosen  angle  in  the  range  0-12'’.  A photograph  of  the  model 
configuration  installed  in  the  tunnel  test  section  is  shown  in  Fig.  3. 

The  flat  test  plate  was  instrumented  with  four  parallel  rows 
of  static  pressure  orifices  located  at  2.92cm  (1.15"),  5.46cm  (2,15"), 
8.64cm  (3.4")  and  11.18cm  (4.4")  from  the  generator  in  the  zero  inci- 
dence position.  Each  of  these  rows  had  4"/  orifices  with  the  spacing 
varying  from  .127cm  (.05")  to  1.27cm  (.5")  depending  on  the  antici- 
pated position  through  the  Interaction.  This  large  number  of  taps 
(i.e.,  188)  was  required  to  ensure  that  complete  detailed  pressure 
distributions  were  obtained  in  the  shock  generator  incidence  range 
0“- I 2“. 

The  surface  pressures  were  sampled  by  a 48  port  Model  48J4 
scanivalve  connected  to  the  model  tappings  via  a bank  of  plastic 
tubes.  The  pressure  levels  were  sensed  using  a DRUCK  strain  gauge 
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type  transducer,  (Model  PDCR22),  digitized  and  then  processed  by 
an  IBM  System  7 on  line  directly  to  an  IBM  158  machine. 
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3.  DISCUSSION  OF  DATA  AND  ANALYSIS 


The  coordinate  system  used  in  presenting  the  data  is  shown 
in  Fig.  4.  Data  taken  along  pressure  tap  rows  parallel  to  the  X 
axis  are  presented  in  terms  of  where  is  measured  rela- 

tive to  the  calculated  shock  location.  This  enables  data  for  dif- 
ferent shock  strengths  to  be  compared  since  it  takes  into  account 
the  movement  of  the  shock  itself  w i t ti  generator  angle.  6|  corre- 
sponds to  the  boundary  layer  thickness  at  the  start  of  the  inter- 
action pressure  rise. 


Boundar\ 


'er  enferinc 


i nteract i on 


As  mentioned  in  the  introduction  the  thickness  of  the  boun- 


dary layer  entering  the  interaction  tends  to  increase  with  increasing 
distance  from  the  shock  generator.  For  the  unswept  leading  edge 
model,  the  variafion  in  boundary  layer  thickness  from  the  inboard 
pressure  tap  row  at  Y = 2.92cm  (1.15")  to  the  outboard  row  at  Y = 

II. 18cm  (4.4")  is  about  66^  at  a shock  generator  incidence  of  2“ 
and  about  44/S  at  an  incidence  of  12°.  The  decrease  in  this  variation 
with  increasing  shock  generator  incidence  is  due  to  the  increased 
angle  of  the  interaction  region. 

In  the  current  test  program  with  the  swept  leading  edge  model, 


a series  of  pitot  surveys  was  made  to  determine  the  boundary  layer 
growth  over  the  new  configuration.  Measurements  were  made  at  several 


streamwise  positions  along  each  of  the  four  rows  of  pressure  tappings. 
Coupling  these  measurements  with  the  pressure  distributions  gives 
the  variation  of  the  boundary  layer  thickness  along  the  start  of  the 
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interaction.  The  variation  in  boundary  layer  thickness  from  the 
inboard  to  outboard  pressure  tap  rows  was  about  50^  at  a generator 
angle  of  2°  and  about  20%  at  an  incidence  of  12°.  This  is  signifi- 
cantly less  than  for  the  original  configuration. 


3 . 2 Surface  pressure  distributions 

Figures  5 through  7 show  the  surface  pressure  distributions 
along  tap  rows  corresponding  to  Y = 2.92cm  (1.15"),  5.46cm  (2.15") 
and  8.64cm  (3.4").  The  pressure  axis  is  the  local  wall  pressure 
non-d  imens  iona I i zed  by  the  freestream  static  pressure.  The  hori- 
zontal axis  is  non-d i mens i ona I i zed  by  6|  , where  6|  is  the 

local  boundary  layer  thickness  immediately  upstream  of  the  start  of 
the  interaction  pressure  rise.  Each  figure  shows  data  for  generator 
angles  from  0°  to  12°  in  2°  increments. 

It  can  be  seen  immediately  from  these  figures  that,  in  terms 
of  6|  , the  upstream  propagation  of  the  disturbance  due  to  the 

shock  wave  is  essentially  independent  of  shock  strength.  Typically 
the  upstream  propagation  is  about  I26|  for  the  two  inboard  tap  rows 
increasing  to  a little  over  I36|  at  Y = 8.64cm  (3.4").  The-.e  conclu- 
sions are  identical  to  those  found  earlier  by  Oskam  (1976)  using  an 
unswept  leading  edge  plate.  As  emphasized  by  Oskam,  this  finding 
is  quite  different  from  two-dimensional  interactions  where  the  up- 
stream effect  is  a strong  function  of  shock  strength.  Depending  on 
the  Y position  of  the  tap  row  and  the  generator  angle,  the  down- 
stream extent  of  the  interaction,  measured  from  the  shock  location 
is  between  20  and  406.  It  can  be  seen  that  1he  larger  interaction 
lengths  are  associated  with  low  generator  angles  and  vice-versa. 
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This  is  brought  out  most^  clearly  for  the  pressure  distributions 
along  the  inboard  tap  row  (Fig.  5). 

The  static  pressure  distributions  obtained  along  the  rows 
Y = constant  have  been  plotted  as  isobar  patterns  and  are  presented 
in  Figs.  8 to  13.  The  pressures  Pj  and  P2  are  the  static  pressures 
in  front  of  and  behind  the  calculated  oblique  shock  respectively. 

Shock  strength  and  shock  angle  have  been  calculated  using  a free- 
stream  Mach  number  of  2.95.  The  pressure  difference  between  Pj  and 
P^  has  been  divided  into  10  equal  increments  and  these  isobars  plotted 
on  the  figures. 

For  a given  generator  angle  it  can  be  seen  that  with  increasing 
distance  from  the  shock  generator  the  isobars  tend  to  spread  away 
from  the  shock.  As  noted  by  Oskam  (1976)  the  pattern  does  not  dis- 
play conical  symmetry  but  does  tend  toward  cylindrical  symmetry 
far  from  the  generator.  Comparisons  with  previous  studies  are 
given  in  Section  4. 

3 . 3 Oil  flow  patterns 

Photographs  of  surface  oil  flow  patterns  for  all  of  the  shock 
generator  angles  are  shown  in  Figs.  14  to  20.  These  were  obtained 
using  a mixture  of  a fluorescent  dye  and  a light  commercial  oil.  A 
very  thin  film  of  this  mixture  was  painted  on  the  surface  of  the 
pressure  plate  id  illuminated  with  an  ultra-violet  lamp.  Once  a 
steady  flow  pattern  had  been  established  at  a particular  generator  I 

angle,  the  photograph  was  taken.  ■ 

i 
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Reference  to  Fig.  14,  which  shows  the  oil  flow  pattern  at  a 
shock  generator  incidence  of  0“  shows  that  the  sweep  of  the  leading 
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edge  does  not  result  in  any  flow  angularity  at  the  surface. 
Throughout  the  entire  region  of  interest  the  surface  streamlines 
are  parallel  to  the  freestream  direction. 

These  surface  flow  patterns  were  reduced  to  quantitative 
data  by  measuring  the  angles  of  local  oil  lines  with  respect  to 
the  X-axis  along  a line  of  Y = constant  (which  coincides  with  the 
pressure  measurements).  Results  obtained  from  that  process  at 
generator  angles  of  2“  through  12“  in  two  degree  increments  are 
shown  in  Figs.  21  and  22.  These  are  at  Y values  of  5.46cm  (2.15") 
and  8.64cm  (3.4").  The  computed  shock  wave  angle  for  each  shock 
generator  angle  is  shown  on  the  left-hand  axis. 

Comparisons  of  the  flow  angles  with  those  of  the  earlier 
study  are  discussed  in  Section  4. 


4.  COMPARISONS  WITH  THE  UNSWEPT  LEADING  EDGE  MODEL 


Figures  23,  24,  and  25  show  comparisons  of  surface  pressure 
distributions  on  the  swept  and  unswept  leading  edge  models  for 
generator  angles  of  4°  and  10°. 

Plotted  in  terms  of  both  the  swept  and  unswept  dis- 

tributions tend  to  collapse  on  one  another  for  any  given  set  of 
conditions.  Sweeping  the  leading  edge  back  by  45°  has  little  effect 
on  the  shape  of  the  pressure  distribution.  In  terms  of  the  thick- 
ness of  the  boundary  layer  entering  the  interaction  in  the  swept 
and  unswept  models  the  difference  naturally  becomes  more  marked 
further  from  the  generator.  For  the  comparisons  along  the  two  rows 
closest  to  the  generator  (i.e..  Figs.  23  and  24)  the  difference  in 
boundary  layer  thickness  is  at  maximum  about  \5%.  Clearly,  dif- 
ferences of  this  order  have  no  significant  effects  on  either  the 
pressure  distribution  shape  or  absolute  pressure  levels. 

For  the  comparisons  at  Y = 8.64cm  (3.4"),  Fig.  25,  the 
entering  boundary  layer  thicknesses  for  the  unswept  model  are 
about  15$  to  20$  higher  at  generator  angles  of  4°  and  10°,  respec- 
tively. Small  differences  in  the  pressure  distributions  can  be 
seen  at  10°.  During  the  initial  pressure  rise  the  swept  leading 
edge  pressures  are  higher  than  the  corresponding  unswept  values. 

Just  ahead  of  the  shock  location,  a cross-over  point  occurs.  The 
trend  Is  reversed  and  becomes  more  noticeable.  In  the  typical 
"dip"  which  occurs  about  56|  downstream  of  the  shock  location, 
the  swept  leading  pressure  are  several  per  cent  below  those  of  the 
unswept  case.  Basically,  the  shapes  of  the  curves  stay  +he  same 


but,  in  the  swept  case,  the  pressure  distribution  is  slightly 
"stretched"  on  the  scale. 

To  more  completely  assess  the  effects  of  the  45®  sweepback 
on  the  whole  interaction  region  (i.e.,  extending  out  from  the  gen- 
erator about  305),  a comparison  of  isobar  patterns  has  been  made. 

These  are  presented  for  shock  generator  angles  of  4®  and  10®  in 
Figs.  26  and  27,  respectively.  For  clarity,  only  those  isobars 
corresponding  to  pressure  rises  of  .2,  .4,  .6,  .8  and  I.O  are  shown. 

Ahead  of  the  shockwave  little  difference  between  the  two  conditions 
is  seen.  The  small  differences  which  are  present,  as  noted  in  the 
pressure  plots,  occur  along  the  outer  tap  rows  where  it  was  antici- 
pated that  the  effect  of  sweepback  would  be  greatest.  Downstream 
of  the  shock  wave  the  differences  in  the  patterns  become  more  marked, 
but  not  iarge. 

To  examine  whether  or  not  a suitable  scaling  para- 

meter the  isobar  patterns  for  the  4®  and  10®  shock  generator  angles 
are  also  shown  in  terms  of  absolute  distance,  X . These  plots  are 
shown  in  Figs.  28  and  29.  Comparison  of  Fig.  26  with  28  and  Fig.  27 
with  29  shows  that  scaling  by  6|  reveals  no  consistent  trends. 

Depending  on  the  generator  angle  and  the  particular  tap  row  examined, 
the  isobar  patterns  may  be  either  drawn  together  or  apart.  For 
several  sets  of  points  along  the  two  outboard  pressure  tap  rows, 
non-d i mens i ona I i z i ng  by  6,  results  in  a measure  of  overcompensation.  | 

I 

It  could  be  that  6|  is  a primary  scaling  parameter  but  that  the  ; 


correct  5 to  use  Is  not  necessarily  the  value  immediately  ahead 
of  a point  in  the  interaction  In  a streamwise  sense.  Since  the  flow 
Is  not  constrained  in  the  Y direction,  as  in  the  two-dimensional 


case,  it  may  not  be  strictly  correct  to  try  and  scale  properties 
along  the  undisturbed  streamwise  direction. 

Comparison  of  the  isobar  patterns  plotted  in  terms  of  X 
(i.e..  Figs.  28  and  29)  with  those  in  terms  of  ^s^^l  (i*®-#  Figs. 

26  and  27)  shows  that  the  spreading  of  the  isobars  with  increasing 
distance  from  the  generator  is  not  noticeably  altered  when  presented 
in  non-dimensional  terms.  This  can  be  seen  most  clearly  for  the  10° 
data.  Both  the  swept  and  unswept  measurements  exhibit  similar  trends. 

Comparisons  of  the  surface  flow  angles  for  the  swept  and  un- 
swept leading  edge  configurations  are  shown  in  Figs.  30  and  31.  Data 
are  shown  for  shock  generator  angles  of  4“  and  10°  at  values  of  Y = 
5.4cm  (2.15”)  and  8.64cm  (3.4"),  I.e.,  about  146  and  246  from  the 
shock  generator.  The  differences  between  the  data  are  small  and 
with  an  accuracy  of  ^1°,  it  can  be  concluded  that  the  45°  sweepback 
has  no  significant  effect  on  the  surface  flow  pattern. 
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Figure  5 Measured  Pressure  Distribution  Y = 2.92cm  (1.15") 


STRTIC  PRESSURE  RPTIO 

.CO  1.25  1.50  1.'5  2.  GO  2.25  2.50  2.^5  2. 00 


i 

I 


ir* 

f' 


‘-iia  DD  -2D.ao  o.m 
XS/DELTP 


20.00  40.00  GO.CQ  SO. 00 

'1=  2.15"  RUN  NO  -21123 


lUG.QG 


FIguro  6 Measured  Pressure  Distribution  Y = b.46cm  (2.15") 


STfiTIC  PRESSURE  RPTIO 

1.50  1.70  2.00  2.2S  2.50 


CJ 


nj 


in 

rj 


rj 

Cj 


li* 

r 

o ^ 


H 


SHOCK  GENERATOR  ANGLE 

.,0  12. 0 


‘-a 


\- 


I 

// 


3-"  -A  1 0 . 0 


\ 


8.0 


G.  0 


rf,  rn' 

^ n 

Tn 

/ i''  P ^- 


>3—  -Eh'- 


.-•ra 


iji'r" 


4. . 0 

2.0 


-f- 


+ 


f,Q.(ju  -HO.  uij  -^n.oc 
yS/DELTPl 


._j.- — 

fiG.GO 


n.cio  20- on  „ Hn.no 

r=  3.4U  RUN  NG  =2424 


80.  OU 
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Figure  8 Surface  Isobar  Patterns 


Figure  11  Surface  Isobar  Potterns 


Figure  12  Surface  Isobar  Patterns  =10 
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